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New possibilities for QCD at finite density 
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I review the growing theoretical indications that at high densities color SU(3) gauge symmetry is spontaneously 
broken by the formation of a quark pair condensate. This leads to a rich phase structure for QCD as a function 
of temperature and chemical potential. 

I also discuss the prospects for lattice QCD calculations at finite density, including the Glasgow algorithm and 
imaginary chemical potential. 



1. Superconducting phases of QCD 

The behavior of matter at high quark density is 
interesting in itself and is relevant to phenomena 
in the early universe, in neutron stars, and in 
heavy- ion collisions. Recent analyses of Nambu- 
Jona-Lasinio (NJL) models of high-density QCD, 
in which the gluons arc replaced by a four- 
fermion interaction, indicate interesting physics. 
In particular, exotic superconducting phases may 
occur above nuclear density. Both two- and three- 
flavor cases have been studied, and I will discuss 
them in turn, before going on to the topic of 
lattice techniques. 



1.1. Two-flavor QCD 

Figure [l] gives a plausible phase diagram 
for QCD with two masslcss flavors. The 
main division is into chirally broken and 
symmetric phases: axial flavor is broken at 
low temperatures/densities, and restored at high 
ones. For a very clear discussion and further 
references see Ref. [[[] . The chiral phase transition 
is believed to be second order at /i = 0, and the 
indications from bag and matrix models arc that 
it is first order at T — 0. Thus there is a tricritical 
point (solid circle) at the switch-over from second 
order to first order [§,[l]. In the low-temperature 
chirally broken region there is the nuclear gas- 
liquid phase transition line, which ends at a 
critical point (empty circle) at T ~ 10 MeV. 

At low temperatures, it is becoming clear 
that additional interesting phases occur above 



the chiral-symetry-restoring chemical potential. 
It was originally suggested by Bailin and Love 
H (see also j§) that QCD at high density 
might behave analogously to a superconductor: 
through the BCS mechanism ||, Cooper pairs 
of quarks condense in an attractive channel, 
breaking the color gauge symmetry, and opening 
a gap at the Fermi surface. Recent mean- 
field/variational analyses of NJL models of QCD, 
using the 4-leg instanton vertex as the effective 
interaction ^,0, indicate that BCS-stylc quark 
pair condensation does indeed occur, and that 
the gaps are phcnomcnologically significant — of 
order 100 MeV — at densities only a few times 
nuclear density. The simplest form of pairing is a 
spin (indices suppressed) and flavor (indices 
singlet, which by antisymmetry of the fermion 
wavef unctions must form a color 3 (indices a, /3): 



(S\qfC^\S) cx 



c a/33 



(1) 



The coherent state \S), consisting of a quark 
pair condensate, has lower free energy than the 
perturbative vacuum, indicating that in the true 
vacuum two quark colors (red and green, say) 
condense, leaving the blue quarks forming a Fermi 
surface. 

In QCD, unlike the NJL model, color is a 
gauged symmetry, so there is no local order 
parameter to distinguish the superconducting 
phase of QCD from the deconfined one. As in the 
standard model, however, we expect there to be a 
first-order phase transition or crossover between 
regions of parameter space with quite different 
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Figure 1. Conjectured phase diagram for QCD 
with two massless flavors. Note the tricritical 
point (black dot), and exotic "superconducting" 
phases at high density and low temperature. 



Figure 2. Conjectured phase diagram for 
QCD with three massless flavors. Unbroken 
global symmetries are given in boxes. The 
superconducting phase breaks baryon number 
and chirality. 



physics. The equation of state, according to 
NJL models, is only slightly affected by quark 
pair condensation. The two main features that 
characterize the superconducting phase are 

(1) There is a gap in the fermion spectrum. 

(2) The condensate changes the electric charges 
of the quarks, since it breaks color and electro- 
magnetism down to color ST/ (2) and a new 
electromagnetism that is a combination of the 
photon and one of the gluons. 

It would be very useful to translate these into 
some observable for heavy ion collisions, but this 
has not yet been done. 

At very low temperatures, a more exotic 
phase may form. The blue quarks, left out 
of the superconducting condensate, may form 
spin-1 pairs and condense, breaking rotational 
invariance and the remaining SU(2) x U(l) em 
gauge symmetry. NJL calculations || indicate 
a parameter-sensitive and small gap (~ 1 keV) 
in this channel, leaving open the intriguing 



possibility that such a phase might play a role 
in neutron stars. 

1.2. Three-flavor QCD 

Fig. shows a conjectured phase diagram for 
QCD with 3 massless flavors ||. We have to use 
the singlc-gluon-cxchangc vertex as the effective 
interaction, since the instanton vertex now has 
an odd number of quark legs, and cannot be 
saturated by a quark pair condensate. However, 
now that the number of flavors and colors is the 
same, a different form of condensate is possible: 

(S\qfC^\S) = Kl 6?8? + k 2 6?6? , (2) 

The mixed Kronecker 8 matrices are invariant 
under correlated vectorial color/flavor rotations 
("color-flavor locking"). The breaking pattern 
is SU(3) color x SU(3) L x SU(3) fl x U(1) B - 
SU(3) diag , where SU(3) diag is the diagonal SU(3) 
subgroup of the first three factors. The color 
symmetry (gauged in QCD) is broken by the 
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quark pair condensate, but, unlike the two-flavor 
case, chiral symmetry is also broken. Also, with 3 
flavors there is a gauge-invariant order parameter, 
corresponding to the breaking of baryon number: 
(NN). 

1.3. Generalities 

The detailed NJL-model calculations that 
show quark pair condensation are given in 
p|. Nambu and Jona-Lasinio explicitly based 
their studies of chiral symmetry breaking on 
BCS, and described the chiral condensate as 
a "'superconductive' solution" ||. However, 
there are important differences (Fig. |^). Chiral 
symmetry breaking is caused by a condensate 
of particle-antiparticle pairs with zero net 
momentum. In the presence of a Fermi surface 
with Fermi momentum pp. one can only create 
particles with p > pp, so as the density grows, 
more and more states are excluded from pairing, 
and chiral symmetry breaking is suppressed. 
In contrast, color symmetry breaking involves 
pairs of particles or pairs of antiparticles. Near 
the Fermi surface these pairs can be created 
at negligible cost in free energy, and so any 
attractive particle-particle interaction enables the 
pairs to lower the free energy. This is the BCS 
instability of the pcrturbative vacuum. If there 
is any channel in which the interaction between 
quarks is attractive, then quark pair condensation 
in that channel will occur. As density increases, 
the phase space available near the Fermi surface 
grows, and more quark pairing occurs. 

There are many directions in which to 
continue investigating quark pair condensates: 
experimental signatures in heavy ion collisions, 
their possible role in neutron star physics, and of 
course the 2+1 flavor case, with a realistic strange 
quark mass. 

Finally, it would be valuable to perform lattice 
calculations of the finite density behavior of both 
the toy models and QCD. Hands and Morrison 
|To| have performed lattice simulations of a Gross- 
Neveu model, which is very similar to the NJL 
model used here, but they have not yet seen 
unambiguous evidence of quark pairing. Finite 
density calculations in QCD, however, remain a 
much more elusive goal. This is the topic of the 
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(a) (b) 

Figure 3. Comparison of chiral and 

color symmetry breaking. Shaded circles are 
particles, open circles are holes (antiparticles). 
(a) Chiral symmetry breaking: particles pair 
with antiparticles of opposite momentum (typical 
pair indicated by arrows), (b) Color symmetry 
breaking: particles pair with particles of opposite 
momentum (typical pair indicated by arrows), 
antiparticles with antiparticles. 



rest of this paper. 

2. Lattice QCD at finite density 

The usual approach to fermions is to integrate 
them out: 

N U(x) sampling weight (3) 

configs 

Sfevm = J x $M?P 

For Monte Carlo evaluation, the sampling 
weight must be positive, so detM must be non- 
negative for any gauge configuration. One way 
to guarantee this is if we have an even number 
of flavors, each with the same fermion matrix M, 
and M is similar to its adjoint, so the eigenval- 
ues are real or in complex-conjugate pairs, and 
det M e R. 

M 1 = PMP^ 1 for some P, N F even (4) 

This is the situation in zero-density lattice QCD. 
For the Wilson action, for example, 

M = 7^ + rD 2 + m + ^ 
Mt = -7^ + rD 2 + m + ^7o 
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so without a chemical potential (|J) is obeyed with 
P = 75 ! but introducing a real chemical potential 
violates this condition, dct M is then complex, 
and straightforward Monte-Carlo methods are in- 
applicable. Similar conclusions obtain for Kogut- 
Susskind quarks. This is the "sign problem", 
which is really a phase problem for QCD. It is 
interesting to note, however, that an imaginary 
chemical potential leaves the measure positive. 

2.1. The Glasgow method 

For the last decade, the approach to finite 
density QCD that has been most seriously 
pursued is the Glasgow method, This 
avoids the positivity problem by treating the 
chemical potential analytically: Z is expanded in 
powers of the fugacity e^, and the coefficients 
are evaluated by Monte-Carlo using the /i = 
weighted ensemble. 



,' det Mfr) \ 
Z(M) -\detM(0)A =0 



(6) 



Unfortunately, the Glasgow method does not re- 
produce the most fundamental property expected 
of QCD at finite chemical potential, namely 
the onset of baryon density when /i reaches 
-Mbaryon/3. Instead, the onset seems to begin at a 
lower chemical potential fi a ~ This is the 

behavior of the quenched theory, in which there 
is an unphysical baryonic pion state |l3( ]. 

In the last year, there has been a convergence of 
opinion that the Glasgow method will eventually 
show the correct onset behavior, but only at very 
large statistics. This is because of the "overlap 
problem" , a mismatch between the regions of 
configuration space that are emphasized by the 
/i = ensemble, and the regions that are relevant 
for the finite-density physics. Because of the 
measure mismatch, the integral comes from large 
but rare fluctuations in the reweighting factor in 
@, making it necessary to accumulate a huge 
number of configurations. 

The evidence for this conclusion has come 
from several sources. (1) Using the Glasgow 
method for QCD on a tiny 2 4 lattice, Barbour 
pl| finds that pL a rises towards mjv/3 at very 
high statistics. (2) The Gross-Neveu model has 
been studied in 2+1 dimensions jL4j. It has 
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Figure 4. Fermion density function of 

chemical potential fi for the Gross-Neveu model 
0, calculated from moderate statistics using 
Glasgow reweighting of ensembles generated at 
chemical potentials //update = 0.0, 0.7, 0.8, as well 
as the exact (unreweighted) result. Note how 
the apparent onset depends on /i up datc, indicating 
that the /z up date = ensemble used for QCD may 
give unreliable results at moderate statistics. 



quarks and conjugate quarks, so there is no sign 
problem, and one can perform updates at any 
chemical potential. At moderate statistics, the 
onset behavior shows dependence on the value of 
fi used in the updating (see Fig. |J). (3) Azcoiti 
ct al |l6| argue that the phase of the quark 
determinant drops as exp(— V), so statistics ~ 
exp(U) is needed. (4) Halasz jlTJ has studied 
matrix models of the QCD fermion determinant, 
and finds that they need statistics ~ exp(iV). 

2.2. Imaginary chemical potential 

As was noted above, with imaginary chemical 
potential fi = iv, det M is real, and the functional 
integral can be evaluated by standard Monte- 
Carlo methods |l8|-|20[|. We still have to choose 
an updating value ^update, and then "reweight" to 
obtain Z(iv), but the reweighting factor is always 
positive, so there is no phase problem, and very 
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high statistics is not needed. 



Ziiv) 



det M(iv) 



Z(wupdatc) \ det M(if updat c) / 



(7) 

with "patches" 



We can cover v = 0...2tt//3 
centered at several different ^update j and thereby 
ensure that the rewcighting factor has arbitrarily 
small fluctuations. Then we must Fourier 
transform to get canonical partition functions 



J N 



2tt 



2tt//3 



(8) 



The Fourier transform is the place where large 
errors may arise. Imaginary chemical potential 
does not bias the ensemble towards states with 
higher quark number; it relies on thermal or 
quantum fluctuations to supply such states, and 
it gives them a characteristic weighting. If such 
fluctuations are too rare then Z(iv) will be dom- 
inated by Zq, and very large statistics will be 
needed to see the effects of the higher Zjy. This 
is not necessarily a problem: as the temperature 
T rises towards the deconfining phase transition 
T c , the baryon becomes lighter, so Mb/T may 
become small enough for thermal fluctuations to 
populate the system with baryons. Calculating 
the Zn by (||) will then be straightforward. The 
Glasgow algorithm may also be expected to work 
better as T -*■ T c . 

2.3. The Hubbard model with imaginary fj, 

As a test of the practicality of imaginary 
chemical potential, it is interesting to explore a 
simple theory that has the sign problem at non- 
zero chemical potential. The natural candidate is 
the Hubbard model in 2 dimensions lllslpoL 



n 



j),<T a i<7 a jO 



(9) 



By a particle-hole transformation, fi — gives 
half-filling. We can replace the four-fermion 
interaction with an auxiliary field A, 

z(n) = ^2 e ~ A2/2 det M 0) det m (-m)- (io) 

A(x) 

The matrices M, given in plf , are real for real /U, 
so the sampling weight is real but not positive, 
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Figure 5. Z(iv) for the Hubbard model on a 
4 2 x 10 lattice, with K = 1, j3 = 1.5, U = 1.0 
pCj ] . The three patches agree well when given the 
same normalization. 



and Monte-Carlo is impossible. For fi zero or 
imaginary, the weight is | det M| 2 , and is positive, 
so we can calculate ratios of partition functions 



Z{iv) 
Z{iv ) 



= ^e- A2 / 2 det|M(^ )| 5 



A(x) 



|detMOo)| 2 ' 



and use several values of ^update = v o to eliminate 
measure-mismatch. In Fig. ^| we show results for 
Z(iv)/Z{ivo) using three values of v . Combining 
them to obtain a single plot of Z(ii/)/Z(0), we see 
that that the three patches agree to within their 
statistical errors. 

Finally, we fit our Z[iv) data to exp(— av 2 ) x 
spline, and Fourier transform it to obtain Zn 
(Fig. ||). At this relatively high temperature we 
have no trouble obtaining Z^ up to N = 5. It 
would be interesting to see whether similar results 
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Figure 6. Zn for the Hubbard model, obtained 
by inverse Fourier transform of the data in Fig. ||. 



could be obtained for QCD at temperatures just 
below the phase transition. 

3. Conclusions 

In summary, while there are growing indica- 
tions from Nambu-Jona-Lasinio and other mod- 
els of exotic phenomena in high-density QCD, it 
remains very difficult to perform the necessary 
lattice calculations. 

The Glasgow algorithm does not appear to 
be practical, since even on a 4 4 lattice it needs 
more statistics than anyone has been able to 
gather. Imaginary chemical potential is an 
interesting alternative, since it is quite possible 
that the numerical problems of inverse Fourier- 
transforming Z[iv) to Z/v will be less difficult. 

Several other approaches are under develop- 
ment, including perfect actions computer 
evaluation of the fermionic Grassman integrals 
j23|, and projecting out the A-particle sectors 
using a new dynamical fermion algorithm p4| . 
Though it is proving hard to find, a successful 
approach will open up new areas of phenomenol- 
ogy to lattice gauge theorists. Clearly the rewards 
will be worth the effort. 
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